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ABSTRACT 


It has been conjectured by Fettweis that wave-digital 
filters designed after doubly terminated ladder networks by 
means of a transmission line transformation have a low 
sensitivity against coefficient variation. The theory of 
wave-digital filters is summarized together with the design 
procedure for deriving wave-digital ladder structures based 
upon an analog design. To investigate this conjecture 
nineteen Chebyshev low-pass filters are realized and 
experimentally evaluated using the IBM-360 generai purpose 
digital computer. The realization was done in two forms: 
aS wave-digital ladder structures and as recursive digital 
filter. The sinusoidal steady state response is determined 
for both filters in floating-point arithmetic with quantized 
coefficients. The analysis shows the validity of the 
conjecture for all nineteen cases over a wordlength range 


of the mantissa of four to twenty two bits. 
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irl N Reber ION 


Transfer properties of analog LC-filters ina ladder 
configuration are known to have exceptionally low sensitivity 
characteristics against element variations. Digital ladder 
structures derived from analog L-C filters retain this low 
coefficient sensitivity (Fettweis conjecture [Ref. 1]). 

This thesis summarizes wave-digital ladder filter theory 
and presents a design procedure for deriving wave-digital 
ladder filter based upon an analog design, and investigates 
the coefficient sensitivity properties of the digital filter 
experimentally. This 1s accomplished by realizing a Chebyshev 
low-pass filter in two forms: aS a wave-digitai ladder 
filter and as a recursive digital filter. The theoretical 
Sinusoidal steady state response for both filter types is 
calculated in floating-point arithmetic with the filter 
coefficients being quantized. The experimental response is 
determined using a computer program on the IBM 360 by 
Carrying out the filter algorithm in floating-point repre- 
sentation under truncation of the mantissa for wordlength 
of four to twenty two bits. A comparison of the results 


demonstrates the validity of the conjecture. 





II. WAVE DIGITAL FILTER THEORY 


eee Ee RINCIPLES 

In conventional digital filter theory filters are 
designed directly from the required transfer function in 
the s-domain so that lengthy analog realization procedures 
can be avoided. Two problems that arise in all digital 
filter design procedures are that of the high sensitivity 
of the transfer function versus coefficients, and roundoff 
noise effects due to finite preciSion arithmetic. Since 
analog LC-ladder structures exhibit low sensitivity for 
element variations Fettweis [Ref. 1] has suggested modeling 
digital filters directly after these analog structures. 
Most of the foilowing introduction is based on that article. 

The basic principle involved is to carry out a frequency 
transtcrmation to convert the system of differential equa- 
tions of the analog structure to a system of difference 
equations that allows implementation by a digital computer. 
The first step in doing this is to perform a Richard's 
transformation [Ref. 2] which converts the analog LC-ladder 
structure to a transmission line circuit by replacing each 
inductor with lengths of short-circuited transmission lines 
and each capacitor with lengths of open-circuited trans- 
mission lines. The length, 2%, of the transmission line is 


related to the inherent delay, by 


Tis 


= i Tr 4 
V 





where v is the propagation velocity. If one assumes a 
traveling sinusoidal signal of wavelength } moving alon 
the transmission line then the period T> between arrival 


of two successive peaks is given by 


bs 


| 
<I> 


By choosing the liength of the transmission line to be equal 


en ee SE ee ee nerd een 


to A one has 


Modeling the transmission line by two parallel pieces of 
delay lines each with an inherent delay t one arrives at 
the unit-element [Fig. 1] which is the basic building block 


for transmission line or also calied microwave filters. 


FIGURE 1. Unit-element of length 2 = i and delay T 


The input impedance to a transmission line of length 2 and 


velocity v is given by 


Z2(s) = Ry tanh (srt) II-4 





wnere Re is the characteristic resistance of the transmission 
line and T = 2&/v. 


Defining a new frequency variable as 
st _ 
Veeee rasa) seein) with 2 27 i= 5 


one arrives at the Richard's transformation. II-5 can also 


be written as 


eel 
a Bes PB Z _ oe ae 
sT Lee Zt 
2 — 2 
for 
Z = aoe IiI-7 
Pee 


Now the input impedance of a short-circuited transmission 


line can be written as 
Z(y) = WR, . II-8a 
and of an open circuited 


os 
ZW II-8b 


10 





In the standard bilinear z-transform, s in T(s) is replaced 


by 





W) 

It 
Fifth 
re 
of- 
|j- 


This is equivalent to the following process starting from 


any analog LC-circuit which has T(s) as its transfer 





mune t Lon : 
fajeebheplace each si. by RK ee ts = R (224) where 
5 a Vv Oo'-z+ 1 
Q 
220 & eS T. eee 
R, = Ll T ane Zee 83 f= 5 
1 SE5 
(b) Replace each sc by R,ctnh (——) 
Q 
eee! 4 ee. mela’ 
where Re - 356 gag z=e ;T= m 


Since resistors are frequency independent elements they 
remain unchanged by this frequency transformation. Thus all 
analog LC-elements are replaced by pieces of open- or 
short-circuited transmission line. 

Fettweis proves very elegantly in Ref. 1 that the 
realization of a digital filter algorithm directly based 
upon the circuit by interconnecting transformed elements 
is physically unrealizable. Also there is no direct way 
of writing an algorithm which takes care of the inter- 
connection of elements. 

Take as an example an inductance. One has to realize 


the voltage-current relation 


al: 





v= VRI . Tee 


Using II-6 to carry out a z-transformation leads to 





v= RI Li—7i 


which results after inverse z-transformation in the following 


difference equation 
v(ntT) + v[(n-1)T] = ae CMU nen) Sere 


Small letters are used to denote instantaneous values. 
It is easily seen that the computation of v(t) at t = nT 
not only requires the knowledge of v(t) at the previous 
sampling instant but also the current at the present one. 
Representing the above difference equation in a signal-flow 
graph in the 7{-domain, one has generated a feedback loop 
without delay which is known to be physically unrealizable. 
This wili be true for the v-1 relationships of all other 
digitized circuit elements. 

This can be avoided by using wave quantities in signal- 
flow graphs and getting from there the algorithm for the 


Gigital filter. For any transmission line the voltage and 


1 





current at any pointcan be defined by a wave traveling 
to the right, A, and a wave, B, traveling to the left after 
reflection. At any instant of time the effective voltage 
is a sum of both waves at any point along the line 

Y= Fetes II-13 
and the current is given by 


me = (A 5) @ tI-14 
O 


Adding both equations gives 


A= > (V + ROI) oes 
and 
p= 2 (v- RI) ET=15b 
2 Oo i 


Leaving off the scaling factor 1/2 gives as the general 
relations between the analog circuit quantities current and 


voltage and the wave quantities A and B 


A= V+ RI II-l6a 


B=V- RI II-16b 
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The next step in the procedure has to be to convert a 
sSignal-flow graph to a wave-flow diagram. By using II-15 
which hoids at any point along the line one calculates 
the v-i relationship which holds only at the input port to 
the transmission line. 


Taking as an example an inductor one has 


V = RI II-17 


with Ro =< we Subst. turIngd Lil santo ll-—)6 and dividing 


II-1l6a and b gives 








pee Vee Yn II-18 
The final step in the process is now to carry out a 


z-transformation of the standard z-transform as given by 


II-7. By substituting 


B= -Z A II-19 
In carrying out an inverse z-transform one finally arrives 


at 


b(nT) = -afl(n-1)T] a 210 
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using small letters to denote instantaneous values. The 
sampling interval is twice the basic delay t. Since a signal 
has to return to the input port after reflection it has gone 
twice through the delay t of the transmission line on a 
round trip to arrive at the next sampling instance. The 


overall process is Summarized in Figure 2. 





RICHARD'S 
TRANSFORMATION Z—TRENSFORM 
ANALOG Se TRANSMISSION a SIGNAL 
LC-LADLER LINE CKT 
VOLTRGE 
INVERSE 
7—TRENSFORM 
7~TRENSFORM ALGORITHM | 
WAVE RELDTIONS UNRFEALIZ2APLE 


/\farcortm 


REALIZABLE 





FIGURE 2. Summary of the steps involved 


Because of the use of waves these digital filters are called 


Wave-digital filters. 
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b. §Gapacitors 


The V-I relation for a capacitor is defined 


as 
R 
V=- I TI~21 
y 


Substituting II-21 into II-16 and dividing II-1l6a and b 


results in 
pe izttazsazia II-22 
The final difference equation is 


bi(nT) = al[(n-1)T] TI-23 





FIGURE 4. Realization of a capacitor 


dee Ciccuie Gdragram 
b. wave-flow diagram 
c. algorithm 
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which leads 


and for the 


maere is no 


what has to 


da. 
by II-16 to 
and 

e. 


Resistors 


Thesv-l relation for a resistor is defined as 


immediately by II-16 to 

B = 0 Milica es 
difference equation 

b(nT) = 0 / bi=Z6 


reflected wave from a resistor if R, = R, and hence 
be required, and hence it is also called a wave sink. 


Open Circuit 


The relation to be realized is I = 0 which leads 
A = B II-27 
a(nT) = b(nT) II~28 


Cherovclrcult 


Pesnout. Cireuit has V = 0 giving by II-16 


A = -B II-29 
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and 


a(nT) = -b(nT) rs 0 


f. Unit-Element 


A unit-element is only defined in transmission 


line filters by the relations 


_ST 
By = Ap = 

_stT BE 3 ayy D 
Bo = Ay e 2 


Note from Figure 1 that the unit-element is used as a two- 


port network. 


The difference equations are 


- _l 
b, (nT) = a, [(n 5) T] 
II-32a,b 


b, (nT) = a, [(n-5)T] 
The delay of T/2 has to be taken into account in determining 
the calculation sequence for the filter algorithm. 
2. Active Elements 
a. Resistive Source 
A source of voltage ve With a series resistance 


R., has a V-I relation 


V. = Vt RI Pi=33 
S Ss 


ane 





where V is the voltage at the output port of this oneport. 


Using II-16 yields 


and 


an). = v,(t) atvti= "nt « 


b. Ideal Voltage Source 


EL 32 


Laie 


An ideal voltage source has the source voltage 


appearing at the output port and hence 


34 


fia 35a5b 


taking R to be an arbitrary constant which will drop out. 


Adding II-35a and b gives 


and 


b(nTtT) = av. (nT) - a(nT) 


20 
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G. £deal Current Source 
Taking the ideal source in combination with a 


series arbitrary R one gets by use of II-16 


B= ae - A II-38 
and 


b(nT) = 2v. (nT) - a(nT) LI-39 


The wave-flow diagrams for the elements of section a. and b. 
are summarized in Table I. 
3. Structural Elements 

So far the different circuit elements have been 
realized. Two problems arise now in interconnecting these 
elements: Reflections from interconnections that are not 
Matched in port impedance and realization of the nodes of 
the analog structure. In the transmission line filter these 
nodes collect and redistribute the forward and reverse 
traveling voltage waves. 

The process of impedance matching is shown best on 


a two-port interconnection 





mn Fa 3 aft + 
0; oy 
V; V5 
= ee ee ee = 
Ry Ro 


FIGURE 5. Interconnection of two ports with 
port resistances Ry and Ro 
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TABLE I 


ELEMENT TRANSMISSION WAV EFLOW WAVEF LOW 
LINE DIAGRAM DIAGRAM 
EQUIVALENT STEADY STATE INSTANTANEOUS 
| L 
|INOUCTOR 


C 


LJ 


CAPACITOR 


R 
——— 
RESISTOR 
UNIT 
ELEMENT 









IDEAL SOURCE 


a 
a 2V. 2V. 
IDEAL SOURC 


Qe ar SS 


HORT CKT 
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At each port II-16 is valid which can be written in closed 


form as 


it 40a, 


Since above shown interconnections are ideal without 
attenuation and delay we have by inspection 
Mei YD 
TIi-4la,b 


PWS TQ, 


Substituting II-41 into II-40 and eliminating tne v's and 


the i's results in 


D5 fg Sy INGE) & in 
pomeecy, * blen “age 
where g iS given by 
Rae i 
Q@ = , II-43 
Ri a Ro 


The wave-flow diagram representation for equation II-42 is 


called an adaptor and the symbol is shown in Figure 6. 


age 








FIGURE 6. Wave-flow diagram symbol of a 


a. parallel adaptor 
b. series adaptor 


Note that due to the fact of Ry > O and Ro > 0 we have 

|8| < 1 and one can see clearly that 8 is the reflection 
coefficient of transmission theory. As shown in Figure 6 
this matching can be extended bo any number of ports. The 
second purpose of an adaptor is to serve as node for the 
wave-flow diagram. 

To arrive at a digital ladder structure of the 
combination of elements one has to generate an algorithm 
that is equivalent to Kirchoff's current or voltase law. 

A series ccnnection of any number of elements has to satisfy 


the relation 


II-44a,b 
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Using IiI-40 for k = 1,n one arrives after use of II-44 and 
elimination cf v's and i's at a system of linear algebraic 
equations that can be written in matrix form. For a three- 


port series adaptor this matrix is 


my P= -8, 1-8, -8. a. 11-45 
bs me “Bz 1-8, ale 


The constant matrix relating output waves to input waves is 


often referred to as scattering matrix. The 8's are given 


by 
2R, | 
68, = = - ITI-46 
k R, +R +... TR, 
and equation in II-45 is of the form 
n 
K k k ee m 


The sum of the coefficients 8 can easily be shown to be 


constant 


@ a Bg te Seo "sp lees 2 TI-48 
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= 


This signifies the presence of only n-1l independent 
coefficients with the nth coefficient as a linear combination 
of all the other ones. 

The corresponding relations for a parallel connection 


of elements are 


Vv, = oma Vay cess Vv TI-49 


eae tesa Je «0 £T=-50 


which upon substitution into II-40 can be shown to lead to 


eueeer Climination of the v's and the i's 


eo er JE LE ey 


and G= 1/R. The sum of the coefficients is again constant 


Pan CLtcteewe se: tt = 2 oie aes, 


aa 2 n 


giving again n-1l independent ones. The matrix equation 


for a three-port parallel adaptor is given as 
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b, a 71 Ot or ay 
bo - Oy Oo-1 a5 ao 1I-54 
b. Oy oD a,—1 as 


eee Ground Rules 
Any realization of a wave-digital ladder filter 
diagram from analog structures has to obey the following 
graund rules: 


1. Connections of ports only pairwise to ensure proper 
wave-flow. 


2. The wave-flow must be compatible, i.e. the outflow 
of one port must flow into the connected port. 


3. The wave flow diagram must not contain delay free loops. 


HO 





III. REALIZATION 


There exist several methods in design of digital ladder 
Structures. A direct deSign was proposed by Bruton [Ref. 3]. 
Another possibility is to synthesize transmission line 
filters by using the theory of unit-element filters [Ref. 4]. 
As a third way one starts from the analog ladder structure 
consisting of reactive elements with resistive termination 
and converts this structure to a digital ladder structure 
by the digitization methods shown earlier. This way was 


chosen in this thesis. 


A. ANALOG DESIGN 

The design of an analog ladder structure was found to 
be quite difficult because driving-point impedance synthesis 
procedures cannot be used to guarantee a maximum power 
transfer over most of the passband.range. The synthesis 
procedure to be used is the insertion loss design procedure 
eege@ee introduction to which can be found in Skwirzynski 
[Ref. 5] where also a great amount of design data is supplied. 
Numerous evaluated filters can be found by Saal [Ref. 6] 
from where the values for the analog ladder filters for 
the experimental realization of the wave-digital ladder 


filter had been obtained. 
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B. DIGITIZATION 

The first step is to replace each analog circuit ele- 
ment in the circuit diagram by its digital counterpart and 
each node by the appropriate adaptor. This procecure leads 
Powcuwmect connections of series - and parallel adaptors which 
generates at the interconnection inner loops without delays 
since any output wave is dependent on any input wave without 
Gelay. A violation of one of the ground rules can be 
avoided if one changes the structure of the analog network 
by use of Kuroda's identities [Ref. 4]. 

With these identities it is possible to convert reactive 
elements from one type to another by introducing unit-elements: 
Unit-elements if matched to the load can be introduced between 
network and load without changing the insertion loss of 
Mme wmetwork (a phase delay will be introduced). By applying 
Rumed@ea Ss icgentities two of which are shown in Figure 8 it 
is possible in a step by step procedure to move these unit- 
elements throucgh the network one by one until all the shunt 


@eeeetrers are converted to series inductors and all inductors 


= 
ny 


are separated by unit-elements as shown in Figures 9, 10 anc i:. 





Figure 8a. One of Kuroda's IGentities 
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Figure 8b. One of Kuroda's Identities. 









Lumped 
Constant 
Ladder 


Figure 9. Introduction of Arbitrary Number 
of Unit-Elements Between Ladder and Load. 


The resistors at both ends in Figure 10 don’t have to be 


equal depending on the particular design. 


IQ 





Figure 10. Application of Kuroda's Identites. 
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Figure ll. Final network with Terminations not 
Shown since they Remain Unchanged. 


This new structure can now be converted to a signal-flow 
diagram without violation of any of the ground rules since 
each adaptor will be separated by a delay represented by 

the unit-element. This realization approach introduces 
unit-elements not contained in the original filter. Since 
they do not generally contribute to the filter action of 

the network their use is not desired. Crochiere [Ref. 7] in 
his work on sensitivity of wave-digital ladder filters uses 
this approach. The main purpose of this section is to 


introduce Kuroda's identities and the use of unit-elements. 


C. TRUE DIGITAL LADDER STRUCTURES 

Another way of solving the problem with the delay free 
loop at the interconnection of adaptors was proposed by 
Fettweis [Ref. 8]. The delay free loop was created due to 
the dependence of any of the output Signals on any of the 
input signals. If one of the two multipliers of the inter- 
Pennected ports is chosen to be equal to 1 then the dcdepen- 


dence of the output signal on the input Signal at that port 


Sal 





is broken and a direct interconnection of adaptors becomes 
possible. No unit-elements are required. This is shown 


on an interconnection of a series - and a parallel 3-port 


adaptor [Fig. 12]. 





Figure 12. Interconnection of a Series 3-port 
and a Parallel 3-port Adaptor. 


The output equation for port R, is given by 


a. + aa. —Jall IlI-1l 


K TI1--2 


and for port Ry by 
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by cela a a,(a, + a. > ac) Pern 


where one has 


2 = III-4 


i eee 


oA 


The delay free locp is easily recognized. Making an arbi- 


trary choice between oe and ys say As one wants 


= Le III-5 


Since R3 is the pert of the parallel adaptor this equation 


is arrived at from III-2 only if 


G, = G, + Gy . c= 


Since the port impedances of ports at interconnection of 
adaptors are not determined by analog circuit element 
values this choice of G3 can be made. The resultant output 


eguation for R3 is 


3 a5 45 + 5 Oo III-/7 


Applying the condition that the sum of all multipliers for 


each adaptor equals two, one can replace a. by (1 - ay) and 


2 


arrives finally at 
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Bee sUMMARY OF STEPS 
(a) By insertion loss design one gets a resistivity 


terminated analog LC-ladder structure [Fig. 15]. 
R3 Rs 
~ vv ——; | 
Vg ii Ro T 3 Te 


Figure 15. Resistivity terminated LC-ladder network 


(b) Apply digitization methods element by element one 


arrives at the wave-digital wave-flow diagram 


PE LG 6.) 








Figure 16. Wave-flow diagram for analog filter of 
Figure 15. 


The port impedances Ry through Ro are predetermined 


by the analog element values. 


(c) Choose the port-impedance for interconnections. 


3) 





Carrying out step three on Figure i6, chooSing in each case 
the right port of each adaptor to have a = 1, (except for 


the last adaptor), leads to the following port impedance 


Gj = G6, + G, III-9a 
Ree =R,, + Ry III-9b 
cy, = G,, + Gy III-9c 
Ras = Ra, + RA III-9d 


The total number of multipliers is six, one for each adaptor 
and two for the last one which corresponds directly to the 
degree of freedom of the analog network (number of resistive 


passive elements minus one). 


E. NUMERICAL IMPLEMENTATION 

To actually implement the wave-digital ladder filter one 
must carefully evaluate the correct computational order to 
ensure the knowledge of ali input quantities for the output 
at hand. Figure 17 shows a conventional signai flow graph 
for the fifth order digital filter of Figure 16 with the 
input and output waves as state variables. For each adaptor 
port two is connected to a delay. From the signal flow graph 
of Figure 17 the computational order can be obtained by 


inspection. Calculations ere assumed to take place 
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Figure 17. Signal flow graph for filter of Figure 16. 
Az) Genotes input to port 1 of adaptor number 3. 


Simultaneously and to take now time. The values at the 

nodes are assumed to be zero between sampling instants except 
for the wave values stored in the delay line. Since an 
actual computer algorithm cannot provide this ideal situation 
one has to follow the signal flow graph. Initial conditions 
‘are assumed to be zero. First the output values for node 
three can be evaluated since it is independent from the 
onmasene input to that port. Going through the signal flow 
graph tc the end one is able to calculate all output values 
for the last adaptor since the resistive termination of the 
last port supplies zero input. Going back through the signal 
flow graph the other outputs can be evaluated. The equations 
for the fifth order filter example in correct computational 


order are: 
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(1) By 3 = 41 (8,478, 9) + ay, 


(VS 5 SEN ey 

Bee a6 “aie a2) * 232 

ee 43 ~ 794, 7 242 

= Bigs = Soils ~ Cabanas 

8) See Say emecey) 

meme si ~ 51°51 * "5.2752 7 gs. 

a 2427 241 5 243 tag lag taytags 

roe 

Bem 4 5 284i) 4 884 88407243? 
Maem 32 8934316 730' 733 

ema = 32490742) © 2327 953 7 237 


LAD boy = Ag, ~ Ang + % (a5, +a50tAaD3) 
(13) boy = ag, 7 plan] +Aagotag3) 

Be a2 ~ 914417412) + 43 

Mea 1242222! * Fn2 F 443 7 F112 


In addition to these equations a storage of output to the 
delays is required. One should also keep in mind the 
equality of waves at directly connected ports, i.e. at any 


instant of time for example b3, is identical to any and az 


1° 


to DD 
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EF. THEORETICAL RESPONSE 
To get the theoretical response it is possible to 
evaluate the transfer function from Figure 17 and let 
Zz = aor A considerably shorter way via matrix theory 
is shown on the fifth order filter. At the interconnections 
of ports in Figure 16 the right hand values of a left adaptor 
are identical tothe lefthand values of a right edaptor. So 
all one has to do is to eliminate port two by algebraic 
procedures which means for a 3-port adaptor: 
(a) Get bo = + za, 
(b) Substitute (a) into the corresponding output 
equation of the adaptor to eliminate bo 
(c) Substitute the resultant aj=in(a,,a) into the 
remaining output equations. 


By changing the dependent quantities in the resultant 


equations of the form 


b, = fn(a, a3) 


1 
IIlfI-lla,;b 
one finally arrives at 
a, = fn(b,,a3) 
III-1l2a,b 


b, = fn(b,,a3) 
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In matrix form, II-12 can be written as 


= III-13 


iL Za 22 


The entries in the P-matrix are generally functions of a 

ame @£ 2. Since b. and a, are identical to the leftside 
values of the next adaptor one immediately generates a chain 
matrix relating ay and by of the first adaptor to b and a3 
of the last one. The last a3 will always be zero (wave sink 
at port 3 of adaptor) so that the chain matrix finally is 


of the form for the example. The matrices for the usual 


=n eee PRI a a || aa ea 
7 OSs 
il Por P22 || 221 “a Roy Roa |} S21 222]! Yo 
circuit arrangements are shown in Table II. Solution of 


the matrix equation will directly lead to the transfer 


b 
ponction H, (2) = = Since the complementary output can 
ial 
be taken from Dia one also is able to evaluate the transfer 


function for the high pass complementary filter 


H, (2) = S 2 
Lal 
; ee tO 
Using z =e one generates the frequency response. 
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TABLE II 







‘a 

CIRCUIT TRANSMISSION | WAVE MATRIX 

SEEMENTS {LINE MATRIX ELEMENTS 
EQUIVALENT 





L\= @ (z+) 
le Pet(zea V/A 
1 Pb Fo (a— | AN 

Be za-1)/A 


& (z-a+il)/A 
i Z| 
Q == (24a 7A 


Q,/° z(a-I)/A\ 
Q =-(Az+ 1I)/A 
Fe 


A= a;(z+1) 
Y « (z- oy ISA 


U = A-z+asI 
2 7S 
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MY. THEGRETICAL SENSITIVITY 


In analog ladder structures resonant frequencies corre- 
Sponding to attenuation poles are determined by the product 
of two element values. The pole and zero positions are 
shifted about their true values with variation of circuit 
element values cauSing a change in magnitude and phase 
characteristic of the response including bandwidth and 
cutoff frequencies. This variation is either caused by 
inaccurate circuit element values or drifting due to various 
physical phenomena. 

In digital filters drifting of values does not exist 
Since the filter coefficients are stored in the computer and 
can be considered constant as long as the filter 1s opera- 
tional. But since the design generates decimal coefficients 
their representation in a binary number of finite wordlength 
generates what is called quantization error. Assuming a 
coefficient obtained with infinite accuracy to be a = 0.82 
its representation in the binary system is (.11010001 ... 

A wordlength of 4 binary digits after the decimal point 
gives 0.1101 which converted to decimal leads to .8125, 
close to .82 but not exact. 

The inaccuracy which is dependent on the number of 
binary digits causes a change in the response. Hence the 
effect of finite wordlength for the digital filter is the 


Same as that of inaccurate analog element vaiues for the 
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analog filter. Analog LC-ladder filters are very insensi- 
tive in their transfer characteristics against element 
variations. It has to be determined whether these sensi- 
tivity properties are preserved and carried over by trans- 
formations. Gupta and Renner [Ref. 9] prove this as follows: 
The voltage-wave into the load in Figure 17 is equal to the 


effective voltage V, over R, since a, = 0 for all times. 





Figure 17. Wave-Digital Two Port 


Since ir e(t) the voltage transfer functions H(s) is 


given by 


i: 
CS) (ei Iv-1 


and the voltage-wave transfer function iS given by 
a 


In(z) | = 2 
HiGz eae 5 eee Tyee 
ay V Ro 


They show the sensitivity of H(s) to be zero against analog 


element variation at particular points, the maximum availabie 
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power points, only. 


SIH(s)i} = 9 Iv-3 
6B 
O 


where fy is the MAP transfer point and §68 is given by IV-10. 


Consider the wave-digital filter multiplier a as a function 


of the analog network elemaents: 


= ty, ey x) Iv-4 


where L,C,R are row vectors. 


A small variation of the nominal value of a will cause 


a corresponding change in L, C and R: 


a+ 6a = £E(L + 6L,C + 6C,R + 6R). Iv-5 


Taking only the first order terms of a Taylor series expansion 


leads to 


© 

Hh 
| 
© 


sn = ane = seh) con) IV-7 


© 
| 
O& 
QD} 
O 
bd | 


or 


6a = K 587 IV-8 
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where 


pec eg gS hy IV-9 
6L 6C OR 
68 = (SLé6CS6R). IV-10 


From IV-l and Iv-2 |H(z)| and |H(s)| are equal: 


|H(z)| = |H(s)]. 


Hence 


6 |H(z) | = 6 |H(s) | Teva 


da 6a 


which upon substitution of IV-8 can be written as 


6 |H(z) | a 6{H(s) | intan 


Sa KR 6e- 


so the zero sensitivity of |H(s)| with respect to variations 

of the analog element values is carried over to the sensi- 
tivity of |H(z)| with respect to variations of provided 

only small variations in the digital-filter coefficients are 
made such that the approximation of the Tayler series expansion 
by the first order term is valid. The zero sensitivity of 


lH(z)| can also be shown in a more qualitative picture 
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Girectly on the wave-digital filter shown in 


Tnhtreoducins the concept of theoretical power (nO Bower is 


Consumed by an algorithm) abscrbed by the wave-twoport one 
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[ete i> < 2 rv-15 


_ ~ -— 1 in ae —- i- + i - . m~m~syes tem 
GCulemime ana no reriections bacK to the source occur, equiva- 


a C4 : ; _~ : = ~ = on «=4 &s ~ =o 
feme tO the MAP points, ome will get the limitinc cése 


‘ EV-Wo . 


THe attenuation A of the wave-digitai filter mev be defined 


as: 





The power P absorbed in the two-port will always be zero. 


So as long as IV-16 is valid at MAP points w_ we have 


O 
A = 0 IV-18 
fase in turn implies 
SA _ g Iv-19 
6a 


One can conclude with thestatement the smaller the attenua- 
tion in the passband the smaller is the sensitivity of the 
wave-digital filter transfer properties with respect to 


coefficient variations. 
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V. RESULTS 


To investigate the sensitivity properties of wave- 
digital filters experimentally 19 different Chebyshev 
low-pass filters were realized. On the basis of analog 
ladder structures taken from Saal [Ref. 7], five filters of 
order N = 5 with ripple (in dB) in the range from R = 1.25 
to R = 0.1 were realized, in addition six filters of order 
N = 7 with ripple in the range from R = 0.28 to R= 0.1. 

The cutoff frequency was normalized to WO, = 1. The sampling 
time was taken to be ie 2.0 for all filters which corres- 
ponds to f/f. = 0.32. The filter with ripple R = .177 

for order N = 5 and order N = 7 was arbitrarily chosen for 

an investigation of dependence of sensitivity properties 

versus sampling time. The sampling times taken were 

T. oer, 1.0, 0.5, 0.1, and 0.05 corresponding to f/f, = 0.32, 
Peo Os, .016 and .008. 

Corresponding to the realization of the wave-digital 
filter a recursive series-cascade digital filter with first 


and second order sections [Fig. 17] was realized directly 


from the transfer function. 





at A nT) 


Figure 17. Series-Cascade Second Order 
Recursive Digital Filter Sections. 
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The choice was made because both filters have the same 
number of multipliers and in general cascade structures are 
known to be less noisy and less sensitive against multiplier 
variations than direct realizations. The design of the 
recursive filter is based upon the transfer function H(s) 
for analog Chebyshev low-pass filters 

H(s) = ————»—— v-1 


1 + e* 1 “(s) 
n 


where z is the n-th Chevyshev polynomial and ¢ is directly 


related to the ripple by 


5 JAS Ve 


eee 120) 


with R denoting the desired ripple in dB and K being a 
constant. Factorization of the denominator polynomial gives 
the pole locations. Collection of complex conjugate poles 


results in the transfer function 


K 


2 2 
(sta) (s +b Stc,)...(s tb ste) 


1 


the relation between pole locations and coefficients in V-3 


is obvious. After application of the bi-linear z-transform 
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one gets H(z) as 





are) 
H(z) = —3— 2 
CZ td, zte,) (z tdoz+te,) +++ (ztd.) 
which can be written as 
Z 2 
si) eae ee zit) : ——Zth) nee 
ro) (Z +d, zte,) CZ, +d,zte.) 
= Ky . Hy (2) cos H (2) IV-5b 


Application of the inverse z-transformation gives the system 
of difference equations that constitute the algorithm of 
the digital filter. 

The response was generated in two ways. After conversion’) 
of the ideal multipliers to binary floating-point representa- 
tion the wordlength of the mantissa was reduced in steps of 
two bits from 22 bits to four bits. The truncated multipliers 
were substituted back into the filter and by use of eo 
the theoretical response was calculated for each step of 
wordlength for both filters, wave-digital and conventional. 
The calculation was carried out in double precision on the 
IBM 360 general purpose digital computer. 

The experimental response was obtained by going through 
Prenaigital filter algorithm up to ten thousand times to 
arrive at the steady-state value. It was assumed that the 


steady state value was arrived at once the normalized output 
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peak had settled down to within 0.001 of the previous peak 
for a sinusoidal input. This was done in steps of Aw = .02- 
within the passband. Multipliers were quantized according 
to the present wordlength of the mantissa and after each 
algebraic step the mantissa of the result was truncated to 
the present wordlength to simulate operation of the digital- 
filter with shorter wordlength for both, memory locations 
and arithmetic element. 

The quantities measured in either case were the maximum 
and the minimum output in the passband, EMGIE)) | eee and 


[H(j) |. , respectively. The experimentally determined 


ripple is given by 


[H(jo) | an 
R ax = 20 log, ———————St‘“s3 V-6 
a |H (3a) | 
max 


The measure of the error in the response is the relative 


error 
Rnax  ®spec 
; = P > 
relative error R Lor oe Rspec 
zZpec 
V-7 
; - < 
relative error 0 EOL Rnax < *®spec 
where ee pec is the specified ripple for the design and denotes 


the maximum variation of the response in the passband. 


ay 





The experimental results confirm Fettweis conjecture 
if one looks just at the theoretical response errors. The 
response after simulation shows no remarkable difference in 
general except for the cases with higher sampling frequency. 
The simulated response is very much affected by the trunca- 
tion noise in the calculation process which is larger for 
the wave-digital filter due to the higher number of additions 


than for the recursive filter, thirty-nine versus fourteen. 
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VI. CONCLUSIONS 


From the observed results it can be concluded that the 
Wave-digital filter is a valuable alternative to conventional 
digital filters. In most cases wave-digital filter can be 
realized with a considerably shorter wordlength than 
conventional digital filters thus allowing an overall 
reduction in size and cost. The larger number of adders 
is insignificant due to the steadily decreasing hardware 
cost and the particular cheapness of adders. The difficult 
access to wave-digital filter due to the unfamiliar insertion 
loss design procedures for the analog LC-ladder filter could 
be eased by generation and publicizing of an insertion loss 
design computer program. Since the sensitivity of wave- 
Gigital filter transfer function properties to coefficient 
Variation does not change noticeably as it is the case with 
the recursive digital filter an application in systems 


with varying or high sampling rates seems very desirable. 
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EXPERIMENTAL RESULTS 


In the following figures the relative error, as specified 
earlier, is plotted versus the number of bits in the word- 
length of filter coefficients and, for the simulation of 
each filter, also in the wordlength of the arithmetic unit. 
Since the range of three decades is not SuEEi cians to display 
all data the results are also shown in tables. An "X" was 
used in all cases where either the filter had become unstable 
or the cutoff frequency was outside the range of observation, 
1.1 radians for a design continuous cutoff frequency of cone 


radian. 
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